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Abstract 

The Coulomb gauge in nonabelian gauge theories is attractive in principle, but 
beset with technical difficulties in perturbation theory. In addition to ordinary 
Feynman integrals, there are, at 2-loop order, Christ-Lee (CL) terms, derived 
either by correctly ordering the operators in the Hamiltonian, or by resolving 
ambiguous Feynman integrals. Renormalization theory depends on the sub- 
graph structure of ordinary Feynman graphs. The CL terms do not have a 
sub-graph sructure. We show how to carry out renormalization in the presence 
of CL terms, by re-expressing these as 'pseudo-Feynman' integrals. We also 
explain how energy divergences cancel. 

1 Introduction 

There are some reasons to be interested in the Coulomb gauge in QCD. It is the 
only explicitly unitary gauge (we discount axial gauges, which have undefined 
denominators l/(n./cQ. It may be useful for bound-state problems. It has 
been used in lattice simulations, and it has been made the basis for a discussion 
of confinement [T], [2], [3], [4], [5], [6], [7], [8], [13], [9], [10], [U (the existence 
of a Hamiltonian allows the use of a variational principle [12]). But there are 
complications in formulating correct Feynman rules, all concerned with the con- 
vergence of integrals over the energy-components of internal energy-momentum 
variables. These problems are eased by using the Hamitonian (phase-space) 
formulation, in which the electric field is one of the dynamical variables 
[T4] . The Lagrangian in the Hamiltonian formalism, and the Feynman rules are 
stated in Appendix A. 

-•^We include here the temporal gauge, in which n is time-like, although this gauge is 
sometimes taken as a starting point. 



Even with this formulation, divergent energy integrals remain in individual 
Feynman graphs, which are cancelled when suitable sets of graphs are combined. 
This has been illustrated in [16 . 

A more subtle complication about the energy integrals remains. This was 
first recognized as a question of operator ordering by [15] (which contains refer- 
ences to earlier related work, and see also [H]). But an equivalent approach is 
to regard it as coming from ambiguous double energy integrals at 2-loop order 
(see [19 ). This integral is (we use bold capital letters for the space components 
of 4- vectors) 

^0 qQ_ 

This double integral is not absolutely convergent, and might be defined as a 
repeated integral in two different ways: 



/ 



^PQ^gO o_^,^ . ■ (1.1) 



p^-P^ ^ze J ^0 - Q 



or 

f dro r dqo -^(ro + go) igp ^ 

The above ambiguous integrals occur in 2-loop order in graphs containing 
just two propagators {0\T{A^{x)Ej{y))\0) (our notation is explained in Ap- 
pendix A) and a connected chain of Coulomb propagators. 

Christ and Lee (and previous authors referenced there) studied the 
correct ordering of operators (including the Fadeev-Popov determinant) in the 
Hamiltonian. They found two new terms called Vi and V2 (of order h^) which 
had to be added. At the same time, the ordinary parts of the Hamiltonian were 
to be Weyl ordered (that is a certain average over different orderings). It was 
argued that this Weyl ordering removed all the ambiguous integrals like (1.1), 
their place being taken by Vi and V2. 

An alternative approach was given in [20^. There, the ambiguity of (1.1) was 
resolved by showing that suitable combinations of Feynman graphs lead to the 
unambiguous, absolutely convergent, combination (defining Po + <7o + ^0 =0) 



r dpo dqo / 
J 2iT 2ix \\ 



ipo iqo , ^qo ^^0 



27r 27r Vpg - p2 + ^2 _ q2 ^ <7o - + ^0 

iro ipo 



r2 - R2 + ie p2 _ p2 + iJ (^-^^ 

This double integral may be done in any order, and therefore confirmed by using 
(1.2) and (1.3). 

Equation (1.4) is related to the identity 

e{ti - t2)e{t2 - h) + e{t2 - h)e{h - h) + e{ts - ti)e{ti - ^2) = -1, 



■^In (section VI) it is stated that the integral (1.7) is zero. There seems to be no 

expHcit statement about integrals like (1.1). If the Weyl ordering implies that all integrals 
like (1.1) are zero, there seems to lead to a contradiction because of the identity (1.4). 
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which has an unambiguous hmit as the times tend to zero. 

Actuahy, (1.4) is not quite sufficient to resolve ah ambguity. In [20 one 
further rule was used. This is that integrals of the form 



f dqodQdrodK , , /(Q)^(R) = (1.5) 



where /, g are any two functions. This rule could not be obtained from within the 
Coulomb gauge, but was deduced from the limit of a gauge which interpolates 
between the Feynman gauge and the Coulomb gauge. The argument is that, in 
the interpolating gauge, the double integrand still factorizes, taking the form 



where each energy integral is now convergent, and so is unambiguously zerojj 
From the identity (1.4) and the rule (1.5), it was shown (provided that the 
external gluons are transverse) in [20 that the additions Vi + V2 of [15] are 
contained within ordinary Feynman integrals. 

The problem with the above treatments of the ambiguous integrals is as 
follows. It is an important property of Feynman graphs that, as a simple ex- 
ample, a 1-loop graph may appear as a sub-graph of a 2-loop graph; and when 
it does it represents the same integral in both cases. This property is vital to 
the process of renormalization of ultra-violet (UV) divergences. The necessary 
counter-terms are found from 1-loop integrals, and are then used to cancel the 
sub-divergences at 2-loops. 

In the presence of the ambiguous integrals (1.1), this property is not obvious. 
For example, the 1-loop graph in Fig. 5(c) contains the integral 



which is naturally, and correctly, taken to be zero. But in the 2-loop graph 
Fig. 6(a), which contains Fig. 5(c) as a sub-graph, the ambiguous integral (1.1) 
appears, which according to [20 has to be combined with other graphs to give 
(1.4). It is not obviously correct to attach the value zero to Fig. 5(c) in isolation. 
What is more, if the sub-graph in Fig. 6(a) is not zero, it is UV divergent, but the 
calculation of Fig. 5(c) provides no counter-term to cancel any sub-divergence 
in Fig. 6(a). In a similar way, the 1-loop graph in Fig. 5(a) is unambiguous (and 
UV divergent), but it appears as a sub-graph in Fig. 6(b) which contains the 
ambiguous double integral (1.1). 

This is a rather trivial example, but in the following we will encounter several 
similar cases. 

It is the purpose of this paper to propose a solution to the dilemma. 
We first have to remind the reader what the functionals Vi and V2 in [15] 
are. 

There is no actual Feynman graph giving an integral like (1.5). There are (non-zero) 
Feynman graphs which give integrals like (1.5) but with ro replaced by po (as in Fig. 8(c)), 
and (1.5) is just something which is added on to these in order to get the combination in (1.4). 




(1.6) 




(1.7) 
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2 The functionals of Christ and Lee 

In order to state what these are, we first need to define the ghost propagator, 
G. It is defined to be the solution of the equation 

{d/dXi) {Df{X)G^''{X, A)) = J^^^^(X - X^. (2.1) 

Here .. are 3- vector indices, a, 6,.. are colour indices, we write a 4- vector 
Xjj, = and Df^ is defined in (A3). In (2.1), A'^{x) is an external gluon 

field which is transverse, dA^{x)/dXi = 0. G is a functional of A^. Note that 
(2.1) is an instantaneous equation, there are no time derivatives. It is convenient 
to define also 

Gf (X, X^ A) = {d/dXi)G^\X, X^ A). (2.2) 

The zeroth order term in the perturbation expansion of G is just the Coulomb 
potential 

G"' = -^"'^ + 0{g) ^ S^'GoiX - X') + 0{g). (2.3) 

In terms of the ghost propagator, the complete Coulomb potential operator 

is 

C"'''(X,X';A) = j d^X"G^i''{X",X;A)Gf{X",X';A). (2.4) 
We need also to define the transverse projection operator 

(X, Y) = SijS^'iX - Y) - didjGo(X - Y). (2.5) 
With this notation, the first Christ-Lee conribution to the Hamiltonian is 

Vi = -(<?V8)r'"r'^''=' J d^Xd^Y 

xGf'iX, Y; A)5^{X - Y)5''''' G^'''{Y,X; A) (2.6) 
For the second, we first define 

Lf^{X,Y;A) = 5''''r^X,Y)+g J cZ^ZGf (X,Z; A)/'=''M^(Z)Tfc,(Z, Y) 

(2.7) 

Then 

V2 = -(578)/"'"=/°''''^' y" dXdYL,"/(X,Y;A)C'''''(X,Y;A)L<'=(Y,X;A). 

(2.8) 

Here we have not used the original form in fTS] but an equivalent given by ^20] 
in his equation (4.5.3). 

The contributions to (2.6) and (2.8) in a perturbation series may be rep- 
resented by 3-dimensonal Feynman graphs, in terms of the Fourier transform 
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(a) 



(b) 



Figure 1: Three-dimensional graphs of order in momentum space, (a) an 
example of a contribution to Vi in (2.6). (b) an example of a contribution to 
V2 in (2.8). The central chain of lines is from C in (2.4) and we call this the 
C-chain. The two side chains are each from L. The zig-zag line denotes (the 
Fourier transform of) Tij. The symbol C on the central segment in (b) draws 
attention to the Coulomb propagator Go in (2.3), which in this case could be 
in any one of three positions on the central vertical C-chain, giving three equal 
contributions. 

A^{K) of A^{X.) (which we assume to satisfy K.A^ = 00- We shah cah these 
CL graphs. Typical examples, to order are shown in Fig. 1(a) for Vi and 
Fig. 1(b) for V2. Our graphical notation is explained in Appendix A. The rules 
are like those for a Euclidean field theory in 3-dimensions. 

It will be useful to have a notation for the CL terms in momentum space, 
like the examples in Fig.l. We will denote these integrals by Vi{N) and V2{N) 
where N is the total number of external gluon lines. Then 



where n stands for the quantum numbers of n external gluons, that is for 
{Kjy^ijy^ajy} wlth u = l,..,n and n -\- n' = N. Thus Fig. 1(a) is an example 
with n = n' = 2. We have put the minus sign on the left in (2.9) because —Vi^2 
are the contributions to the effective action.) 

^This condition is assumed also in the work of Christ and Lee. Thus, unhke the contribu- 
tions form ordinary Feynman diagrams, the CL part of the effective action is not known for 
non-transverse external gluon fields. 





(2.9) 



5 



Similarly, 

-VoiN) = - fabcra b c 

= ^ I dWCiW2{P,Ci,'R), (2.10) 

where here n + + n'^ = TV. Fig. 1(b) is an example with n = = 1, n^' = 2. 
The factor (n^^ + 1) is inserted into (2.10) to take account of the n^' + 1 possible 
positions of the zero-order Coulomb propagator (2.3) within C (as indicated by 
the 'C in Fig. 1(b).) In (2.10), P is defined to be the momentum through the 
transverse projector {Sij — (P^Pj/P^)} in L(P), and Q likewise. R is defined 
so that P + Q + R = 0. 

The theorem proved in [20 may be stated as follows: the sum of the am- 
biguous Feynman integrals (containing (1.1)) may be brought to the form (after 
some manipulation) 



d'^p d^q / Po ^0 ^0 Po_ 

(27r)4 (27r)4 Vpg - ^ ie - ^ ie - ^ ie p^ - P^ 



^0 Qo 



R2 + g2 _ Q2 + 



X (2.11) 



where W^^ and are defined in (2.9) and (2.10) and X is an integral of 
the form of (1.5), and which is taken to be zero. The integrals in (2.11) are 
4-dimensional; but not all the terms are Feynman integrals generated from the 
rules in the appendix, because the denominator in (1.7) can appear without the 
transverse projector {6ij — (PjPj/P^)} in the numerator. Also, since (1.4) is 
independent of P, Q, R, the identification of these with the integation variables 
in (2.9) and (2.10) is arbitrary. 

If dimensional regularization is used, with 3 — e spatial dimensions, the 1- 
loop sub-integrations in CL graphs like those in Fig.l have no divergences, that 
is no poles at e = 0. The complete 2-loop integrals are UV divergent, with 
single poles at e = 0. At first sight, it is unexpected that there are no divergent 
sub-graphs in CL graphs. But these sub-graphs have no meaning on there own; 
there are no 1-loop CL graphs. 

As explained in section 1, the problem is how to combine the CL graphs 
with ordinary Feynman graphs, particularly for the purpose of renormalization. 
In the next section, we propose an answer to this dilemma. 

^This problem is briefly alluded to in 15 at te end of section VI. It is noted that calculation 
of an order g"^ contribution to Vi involves singular terms of the form (5^(X)/|X|. It is suggested 
that these may be 'relevant to the cancellation of divergences from the usual two-loop Feynman 
graphs'. Using dimensional regularization, we do not encounter such singular terms. 
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3 Expressing Christ-Lee terms as Feynman-like 
integrals 



We want to make (2.11) more like ordinary Feynman integrals. For the 
term this is easy, but for the Wi' term it can only be achieved partially. 

The order of the po^Qo integrations in (2.11) is optional, because of the 
identity (1.4). Let us choose the order to be that in (1.3). Then only the 
first term in the square bracket (2.11) contributes; the other two terms are 
zero because of equations like (1.2). Then, in the contribution from VKj^, there 
appears the correct Feynman propagtor 



and similarly for Q. An example is shown in Fig. 2 (a). This looks like an ordinary 
Feynman graph, but it is supplemented with the rule that the ambiguous integral 
(1.1) is in this case to be interpreted in the order in (1.3). Thus we are to some 
extent reversing the process carried out in [20] . But now we attribute the value 
of (1.1) to a single type of graph, rather than to the combination in (2.11). 
The instantaneous propagators (denoted by dotted lines) are the same in the 
graphs in Fig. 2 as in the CL graph in Fig. 1(b). We will call graphs like Fig. 2(a) 
'pseudo-Feynman graphs'. 

The above process does not generate any graph like the one in Fig. 2(b), 
although this looks like an ambiguous Feynman graph. Thus we may say that 
pseudo-Feynman graphs like Fig. 2(b) are zero. 

Note also that the above rule is consistent with the term X in (2.11) being 
zero (because of (1.2)). 

The general rule is that the energy integrals are to be done in the order with 
the energy on the C-chain of lines (defined to be a factor like C as defined in 
(2.4)) last. The C-chain is the central one in Fig. 2(a) but the right hand chain 
in Fig. 2(b) (it is the chain which contains a segment with just the Coulomb 
propagator 1/K^)). However, there are some special exceptions to this second 
rule, arising form Vi. 

We have not been able to express all of the Wf in (2.11) in terms of pseudo- 
Feynman graphs, but just some special parts of W^^ . But these will turn out 
to be sufficient to treat renormalization. Fig. 3(a) shows an example of a class 
of diagrams of the form of the CL graphs in Fig. 1(a), which have any number 
of external gluon lines on the left, but only one on the right. The CL integral 
has the form 



where the details of the function F are not important (and n^' is the number 
of external gluons in F). Because the external gluon on the right is transverse, 
we can replace Q[ by —Qi. We can also insert a term Sik J d^Q{l/Q''^) which 




(3.1) 




d^Q Qk Q'i 
(27r)3 Q2 Q'2 



(3.2) 
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Figure 2: (a) is a pseudo-Feynman graph which is equivalent to the CL graph 
in Fig. 1(b), provided that the energy integrals in the inserted factors are done 
in the order in (1.3). (b) is a Feynman graph of a similar type, but the C-chain 
(which in this case reduces to just a single line) is now the one on the top right, 
and the energy integrals are to be done in the order of (1.2), consequently giving 
a zero contribution. In this graph, r' = r + /ci so that p + g + r = 0. The symbol 
C draws attention to the propagator Go in (2.3). In (a), this segment can be in 
any of three places, giving three equal contributions. 



is zero with dimensional regularization. Thus we get instead of (3.2) 

Now we insert the integrand in (3.3) into (2.11), and obtain 

(2,)- / A-MK,^_ - ^) (3^4, 

where the po and ro integrals are to be done in that order. (3.4) is almost a 
Feynman integral, except that the p-propagator has only the numerator 5ik not 
the transverse projection operator in (3.1). This incomplete Feynman integral 
is denoted by the graph in Fig. 3(b), where the cross bar on the p-line denotes 
its incomplete numerator. Apart from this feature. Fig. 3(b) is a particular case 
of the class in Fig,2(b). Thus there are special exceptions to the statement that 
the terms in Fig. 2(b) are zero. 

(Colour factors are omitted in (3.3) and (3.4).) 

4 Renormalization 

Renormalization counter-terms are derived from 1-loop integrals. These have 
to cancel divergences in sub-graphs at 2-loops. The problem is how to identify 
these sub-graphs when there are CL terms which do not have the sub-graph 
structure of normal Feynman graphs. We will show that, re-expressing the CL 
integrals as pseudo-Feynman graphs, like those in Figs. 2 (a) and 3(b) enable us 
to solve this problem. 
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q 

(a) 



.y 

(b) 



Figure 3: (a) A particular class of CL diagrams from Fig. 1(a), with only one 
external gluon on the right, (b) An incomplete pesudo-Feynman diagram which 
gives something equivalent to (a). The cross bar on the p-line denotes the 
incomplete numerator. 



Q O 




Figure 4: 0{g^) graphs containing 0{g'^) counter-term insertions, which are 
represented by small circles. 



The divergent part of 1-loop effective action is known, for example in pT] : 



-Cg 



-2 1 r II 



167r2 



Oi 



{FSddjAI - gr'^A'iA'^^} + (l/2)Ef + Efd^A^o + K + diC*")diC^) 

(4.1) 

where F^jiv is defined in (A2), Cq is the Casimir of the colour group, c is the 
ghost field, and the source is introduced in (Al). The Feynman graphs and 
pseudo-Feynman graphs together have to have the right UV divergences to be 
cancelled by these counter-terms. In Fig. 4 we show some 0{g^) graphs with 
counter-term insertions deduced from (4.1). For our present purposes, the only 
relevant term in (4.1) is the second one. 



(a) 



Figure 5: 1-loop self-energy graphs, (a) and (b) are Coulomb line self-energy 
graphs. 



P 




Figure 6: 0{g^) graphs containing self-energy subgraphs, (c) is an ordinary, 
unambiguous, Feynman graph, (b) is a pseudo- Feynman graph of the type of 
Fig. 2(a), and Fig. 6(a) is of the type of Fig. 3(b). According to the prescription, 
in Fig. 6(a) the integrals are to be done as in equation (1.5), giving zero; but 
in Fig. 6(b) the integrals are to be done as in equation (1.3), that is integrating 
over the self-energy sub-graph first. 

We begin with Coulomb self-energy sub-graphs. The 1-loop graphs which 
give the counter-terms are shown in Fig. 5 (a) and (b). There is an important 
cancellation between these two graphs, which ensures that their sum is propor- 
tional to K^, which is necessary in order to cancel one of the two denominators 
from Coulomb propagators. Fig. 6(c) is a 2-loop Feynman graph which contains 
Fig. 5(b) as a sub-graph. Fig. 6(b) is a pseudo- Feynman graph (equivalent to 
a contribution from the V2 CL term) which contains Fig. 5 (a) as a sub-graph. 
Graph 6(b) comes with the prescription that the po-i^tegral is to be done first, 
then the ro-integral; and this is just what is required for renormalization, where 
divergent sub-integrations are to be done first. 

Next we turn to UV divergent triangle sub-graphs. Fig. 7 shows two such 
graphs. Fig. 7(b) is an ordinary, unambiguous Feynman graph, in which p, 
and q lines make up a UV divergent triangle sub-graph. Fig. 7(a) is a pseudo- 
Feynman graph, of the class of Fig. 2(a), also with a UV divergent sub-graph 
triangle. In Fig. 7(a), the po-integral is to be done before the ro-integral, and 
this is the order needed for for the cancellation of the UV divergence (by the 
counter-term in Fig.4(e)). 

Fig. 8 shows three more graphs, in which the triangle formed by the p, q' 
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Figure 7: 0{g^) graphs containing UV divergent triangle sub-graphs, made up 
of the hues p^p'^q. (b) is an ordinary, unambiguous Feynman graph, (a) is a 
pseudo-Feynman graph, of the class in Fig. 2(a), derived from a CL term. In 
(a), the po-integral is to be done before the tq one. 






Figure 8: 0{g^) graphs containing UV divergent triangle sub-graphs made up 
of the lines p^q^q' . (b) is an ordinary, unambiguous Feynman graph, (a) is 
of the class in Fig. 2(b) and so is zero, (c) is derived from Wi in (2.11), and 
cancels the UV divergence in (b). In each graph, the symbol C draws attention 
to the position of the Gq propagator in (2.3). The bar on the p-propagator in 
(c) indicated that only the 5ij term in the numerator is retained. 



lines is an UV divergent sub-graph. Again Fig. 8(b) is an ordinary unambiguous 
Feynman graph. Fig. 8(a) is of the class in Fig. 2(b) and so gets zero contribu- 
tion from W2 in (2.11). But there is a relevant graph Fig. 8(c) of the class of 
Fig. 3(b), coming from Wi in (2.11). The difference between (a) and (c) lies 
in a factor PiPj/P'^ and because of the transversality of the q propagator this 
can be replaced by PiR^/P'^, and this contributes only an UV finite part to 
the sub-integration. Thus the UV divergent sub-integrations in (b) and (c) do 
combine as they should, and in fact cancel each other; so there is no need for a 
counter-term. 

We have verified (at 0(^^)) the existence of sub-graphs with UV divergences 
necessary to cancel the counter-terms in Fig. 4(a), (b), (c) and (e). In the cases 
of (d) and (f), no ambiguous graphs are involved. 
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5 Energy divergences 



Another unfortunate feature of the Coulomb gauge is the presence, in individual 
graphs, of energy-divergences, that is integrals of the form 

J droF(ro,R) (5.1) 

where 

F F(R) as ro ^ oo. (5.2) 

Some, but not all, of these energy divergences occur in just the same graphs 
which have the ambiguous integrals described in section 1, that is (to 2-loop or- 
der) graphs containing only two transverse gluon lines. The ambiguous integrals 
show up when doing the integrals in the order 

J d^Kd^Pdrodpo. (5.2) 

The energy divergences appear from the order 

J d^Kdrod^Pdpo (5.3) 

which is the order required by renormalization theory. 

It is expected that these energy divergences should cancel when sets of graphs 
are combined. This has been confirmed at 0{g^) in [17] but only when the sub- 
graphs are quark loops. Then, the 0{g'^) effective action is gauge- invariant, 
obeying Ward identities. The analysis in [TT^ made use of these Ward identities. 
We want to extend the argument to cover gluon loop sub-graphs. Then the 
0{g'^) effective action is in general not gauge- invariant, but only BRST invariant, 
so it is not obvious that the cancellation of energy divergences is as simple. The 
divergent parts are shown in equation (4.1), and indeed they are not all gauge 
invariant. However, the only one relevant to the energy divergences is the second 
term in (4.1), which is gauge invariant. 

The counter-term graphs in Fig. 4 exhibit the energy divergences clearly, and 
the 2-loop graphs whose sub-divergences Fig. 4 cancel are the ones with energy 
divergences. 

In order to extend the argument of [17] to gluon loop sub-graphs, we need 
to show that these sub-graphs at high energy obey Ward identities, not just 
BRST identities. The divergent parts obey Ward identities because the ghost 
graphs shown in Fig. 9 are UV convergent. In order to extend the argument 
to high-energy limits, we need to verify that the graphs like those Fig. 9 are 
also suppressed at high energy. The graphs in Fig. 9 have open ghost lines 
terminating at the ghost field and the sources Uq^u^ or ^ which appear in 
the original Lagrangian in equation (Al) in Appendix. 

The graphs in Fig. 9 all appear by power counting to be linearly divergent 
in the UV, but they are each proportional to Pi (because the gluon line k is 
transverse), and so are at most logarithmically divergent. But in (b) there are 
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Figure 9: Graphs contributing to the BRST identities. These graphs contain 
open ghost hues, terminating at the Zinn- Justin source Uq^u^ or ^ and the 
bottom hue is a ghost hue. The integrals appear by power counting to be 
hnearly divergent in the UV, but they ah develop external spatial momenta in 
their numerators which cause them to be UV convergent, (b), (c) and (d) are 
energy-suppressed but the example in (a) is not. 

three spatial indices, /c and so there must be another external spatial mo- 
mentum. And there is further cancellation between (c) and (d), which produces 
an extra power of external spatial momenta. Thus (b) and (c) + (d) are actually 
UV convergent (as indicated by the absence of trilinear ghost counter-terms in 
(4.1)). Also, to balance dimensions, we expect at least one power of po in the 
denominator, suppressing the integrals at high-energy (that ispo ^ IPI, IQI? <7o)- 
In Appendix B, we show this suppression in more detail. 

The above suppression does not occur in all ghost graphs. Fig. 9(a) is UV 
convergent, but proportional only to one power of external spatial momenta 
(Qi), and is in fact independent of and q^. But, unlike (b), (c) and (d), graphs 
like (a) with or Uq sources do not enter into the BRST identities multiplied 
by a factor of energy, and so are irrelevant to the high-energy behaviour of the 
vertex parts. 

6 Conclusion 

To two loop order, the existence of the Christ-Lee terms in the Hamiltonian 
for Coulomb gauge QCD seems to be a problem for renormalization, because 
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they do not have sub-graph structure of ordinary Feynman graphs. We have 
shown that this difficulty is overcome, by re-expressing the Christ-Lee terms as 
pseudo-Feynamn graphs. We have also shown how energy divergences cancel 
between Feynman and pseudo- Feynman graphs. 

However, to three-loop order there may remain problems about the renor- 
malization of the Christ-Lee terms, when ultra-violet divergent sub-graphs are 
inserted into the ambiguous 2-loop graphs. This question has been discussed in 
section 5(ii) of |20j and in fST, [22 . 

This paper is only about perturbation theory. An important open question 
is whether there are related complications in, for example, lattice calculations. 
Certainly, the Hamiltonian must be correctly ordered, and this generates CL 
terms. 



Appendix A - Feynman rules 

The Lagrangian density, in phase space form, which we use is 

where fi = {o^i)^ c is the ghost field, u^^u^^vf and are sources. For a 7^ 0, 
the equations of motion are non-singular, and we can derive Feynman rules. 
Then the limit a ^ is to be taken, to get the Feynman rules for the Coulomb 
gauge. In (Al) 

and 

Df = S'^'d^ - gr'^A^ (AS) 

In (Al), Aq appears linearly and with no time derivative, so it may be 
integrated out; although we state Feynman rules in the form in which Aq is 
retained. If it is integrated out, it enforces the Coulomb law 

DfE^^ = 0. (A4) 

From (2.1), this gives 

EtiX) = Ef^iX) + J d^YGfiX, Y; A)p''(Y) (A5) 
where diEf"- = and 

= gr'^AlEj-. {A6) 

Using these equations and (2.4), 

^ jd^XEfEf = ^ j d^XEl'^El'' + ^ J d3Xd3xV°(X)C°''(X,X';A)p''(X') 

(AT) 
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Figure 10: The Feynman rules for propagators. The general rule is that dashed 
lines represent A^, dotted lines represent and continuous lines represent E*^ . 



The second term in (A7) is the non-abelian Coulomb potential operator. 
How to order the operators (which are all at the same time) in (A7) is one of 
the questions which lead to the CL terms Vi and ¥2- 

The Feynman rules for (Al) are derived for instance in [17 . They are shown 
in Fig. 10 and Fig. 11. In the Coulomb gauge, ghost lines are identical with parts 
of Coulomb lines. In fact closed ghost loops just cancel closed Coulomb loops. 
Therefore, only open ghost lines are required, where the ghost is attached to 
one of the sources u^^Uq^v'^ ^ which appear in (Al), and are shown in Fig. 11. 



Appendix B - High energy limits 

Here we examine the high energy limit of ghost Feynman graphs in more detail 
than in Section 5. We take as an example Fig. 9 (c) and (d). The sum of these 
two graphs contains the integral 

P f d^KdJz,, ^'^^^^ TiijK') 1 2 J, . . 
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Figure 11: Feynman rules for vertices. The last four involve open ghost 
(with colour index c) and ghost sources. 
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Doing the ko integral gives 

Hqo ^ -go). (^2) 

We are concerned with the hmit qo |P|, |Q| and this comes from the region 
of integration |K| |P|, |Q|. So the high-energy Hmit of (B2) is 

IT / ^'^(K^ ( '^]k|'^' ) ^^■'(K)^'^(K') + ^ + ''(1/90). 

(S3) 

provided that this K integral is convergent. It is convergent because the factor 
(|K'| — |K|) provides at least one power of Q and by rotational invariance there 
must be another factor of P or Q. So (Bl) tends to zero faster than I/qq at 
high energy. 

By a similar argument, the graph in Fig. 9(b) is shown to have a high-energy 
limit of the form 

{l/qo)QjMju{P,Q) {B4) 
where M is a third rank tensor with dimensions of a momentum.. 



We are grateful to Dr. G. Duplancic for e-drawing the figures. 
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